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Strain-induced modulation of magnetic interactions in graphene 
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The ease with which the physical properties of graphene can be tuned suggests a wide range of 
possible applications. Recently, strain engineering of these properties has been of particular inter- 
est. Possible spintronic applications of magnetically-doped graphene systems have motivated recent 
theoretical investigations of the so-called Ruderman-Kittel-Kasuya-Yosida interaction between lo- 
calized moments in graphene. In this work a combination of analytic and numerical techniques are 
used to examine the effects of uniaxial strain on such an interaction. A range of interesting fea- 
tures are uncovered depending on the separation and strain directions. Amplification, suppression 
and oscillatory behaviour are reported as a function of the strain and mathematically transparent 
expressions predicting these features are derived. Since a wide range of effects, including overall 
moment formation and magnetotransport response, are underpinned by such interactions we predict 
that the ability to manipulate the coupling by applying strain may lead to interesting spintronic 
applications. 

PACS numbers: 



Graphene has been attracting the attention of 
the wider scientific community due to an enormous 
range of tuneable properties, suggesting applications in 
fields as diverse as photonics, sensor technology and 
spintronics tlna ■ In recent yea rs, the p otential to tune the 
electronic [ilnLtransport [i,[lM3|, optical [l3,[ll,[i3| 
and magnetic |2Cll - l22| properties of graphene systems by 
applying strain has been explored. The degree to which 
these properties can be tuned is enhanced by the different 
types of strain that can be applied. Apart from simple 
uniaxial strains 0, more exotic features like creases 
and bubbles can be introduced 0, [ill, [H, 0, IH IHl . 

An important topic in spintronics is the indirect 
exchange interaction between localized magnetic mo- 
ments mediated by the conduction electrons of a con- 
ducting host medium. This interaction manifests it- 
self as an energy difference between different align- 
ments of the localized moments, leading to energeti- 
cally favourable configurations. Such an interaction is 
usually calculated within the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) approximation j25 - 2 7 and indeed the in- 
teraction itself often takes this moniker. The RKKY in- 
teraction in graphene has been intensively studied [28- 
I43I I with a general consensus that the interaction streiigth 
decays asymptotically as D^^ in undoped graphene |44l|. 
where D is the separation between the magnetic mo- 
ments. This fast decay rate, arising from the graphene 
electronic structure at the Fermi energy, results in the 
interaction being very short-ranged. Any method of am- 
plifying the coupling to extend its range is welcome and 
could prove useful for future spintronic applications. An- 
other peculiar feature of this interaction in graphene- 
based materials is the masking of the usual sign-changing 
oscillations due to a commensurability effect [45i] . 

With the motivation of amplifying the magnetic inter- 
action strength, in this letter we show how the magnetic 



coupling between localized moments in a graphene sheet 
can be manipulated by applying uniaxial strain. The 
sequence adopted is as follows. We start by introducing 
the general formalism used to calculate the magnetic cou- 
pling, which is written entirely in terms of the real-space 
single-particle Green functions (GF) of the host graphene 
sheet. We subsequently show in a mathematically trans- 
parent form how the GF is affected by the applied strain 
and use this result to predict the behaviour of the cou- 
pling when the direction or strength of the strain is var- 
ied. We find that both amplification and suppression 
of the magnetic couping can be achieved. Furthermore 
we demonstrate that inter- and intra-sublattice couplings 
can be switched on and off independently, suggesting a 
wide range of possible applications. Our results are then 
confirmed using a fully numerical calculations. 

Methods We start by considering two substitutional 
magnetic impurity atoms at sites A and B a distance D 
apart embedded in a graphene sheet. Despite the sim- 
plicity of this setup, it is sufficient to capture the essence 
of the magnetic interaction between arbitrary magnetic 
objects. The indirect exchange coupling between these 
two moments can be calculated by considering the en- 
ergy difference between the ferromagnetic (FM) and an- 
tiferromagnetic (AFM) alignments of the moments (46| . 
This energy difference can be written as a perturbation 
expansion in powers of the exchange splitting Vex be- 
tween the spin-polarized orbitals of the magnetic impu- 
rities and when expressed to leading order in Vex gives 
an expression equivalent to the commonly used RKKY 
approximation 
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dEf{E)Im[gBA{E)gAB{E)], (1) 



where gBA{E) is the single-electron GF describing elec- 
tron propagation in the pristine host material and f{E) 
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FIG. 1: Schematic representation of the graphene lattice with 
the armchair (zigzag) direction marked by the arrow labelled 
'A' ('Z'). The filled and hollow symbols represent sites on dif- 
ferent sublattices. The unit vectors ai and 3.2 and unit cell 
(dashed rectangle) are also shown. The large (red) symbols 
represent magnetic impurities in the lattice, in this case sep- 
arated by the unit of separation in the armchair direction 
D = Ia- The unit of separation in the zigzag direction, Iz, is 
also shown. The distances R\, R2 and -R3 between an atom 
in the lattice and its three nearest neighbours are shown. 



is the Fermi function. We note that the moment separa- 
tion, £), enters only in the GF terms. 

To calculate the GF of the host graphene sheet, 
a nearest-neighbour tight-binding Hamiltonian is used. 
The electronic structure of graphene is well described 
using this model and furthermore it has been ex- 
tended successfully to deal with strained graphene 
systems Q- A Hamiltonian of the form % = 

Sr,r' Sn,n' C,'^ IS uscd, where |r,n) labels 

a TT orbital centred at site n = or 1 in the two-atom unit 
cell shown in Fig. [T] whose location is given by r, t^'^, is 
the electronic hopping term between two such orbitals. 
The summations are over nearest neighbour sites only. 
Each of the orbitals considered has three such neighbours, 
again shown in Fig. [TJ In unstrained graphene, the bond 
lengths i?i, i?2 and i?3 are identical, and therefore so are 
the associated hopping terms ti , ^2 and ^3 which take the 
value io = — 2.7eV. When a tensile strain, e, is applied 
to the graphene sheet the bond lengths and hence the 
hopping values are altered. 

For uniaxial strains in the high symmetry zigzag (Z) 
and armchair (A) directions the bond lengths vary with 
strain as follows : 
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bond length in graphene and a = 0.165 is the graphite 
value for Poisson's ratio, giving the level of contraction 
in the direction perpendicular to the applied strain. The 
hopping parameters vary with the change in bond length, 
AR as 



t{AR) = to e" 



(3) 



where a = 3.37 is taken from the literature [1, |43|. Us- 
ing the hopping parameters found using Eqs. ^ and 
([3]) allows us to calculate the bandstructure of a strained 
graphene system. The dispersion relation is given 
by e± = ±1/^2 + 4 tit2 cos kA cos kz + 4:tl cos^ kz ■ For 
convenience we have defined dimensionless A:-space vec- 
tors kA = ^ Ia kx , kz = ■^Iz ky in terms of Ia Hz) - 
the strained unit of length between unit cells separated 
in the armchair (zigzag) direction and shown in Fig. [TJ It 
is important to distinguish between the strain and sepa- 
ration directions. Both armchair and zigzag separations 
between the moments will be considered, and for both 
cases strains will be applied parallel and perpendicular 
to the separation direction. In our convention, zigzag is 
perpendicular to armchair, so that the parallel (perpen- 
dicular) strain direction corresponds the same (opposite) 
high symmetry direction to the separation. 

The real-space GF between two sites on the graphene 
lattice separated by a vector D can be written as a dou- 
ble integral over the Brillouin Zone. We have shown 
previously [3^ for unstrained graphene that one of the in- 
tegrals can be solved analytically using contour integra- 
tion and that for high-symmetry direction separations, 
the remaining integral is very well approximated using 
the Stationary Phase Approximation (SPA). This ap- 
proach allows us to write the GF for energies throughout 
the entire band in the form 



Gd{E) 



A{E)e 



iQ(E)D 



D 



(4) 



where A{E) is an energy-dependent coefhcient and Q{E) 
can be identified with the Fermi wavevector in the direc- 
tion of separation. The exact functional forms of these 
quantities depend on the separation direction, but the 
distance dependence of the GF is clear in this form. Fol- 
lowing Ref. [3^ we can generalise the expressions to 
strained graphene. 

For separations in the armchair direction between sites 
on the same sublattice, and for energy values in a broad 
range around E = Q {\E\ < 0.5|to|), we can write 



AiE,s)^J- 



QiE^e) 




(5) 



where, i?i = R3 due to symmetry, Rq is the unstrained 



For zigzag separations there are two contributions (±) to 
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FIG. 2: Comparison between numerically calculated (sym- 
bols) and analytic expressions (lines) for the GF between two 
sites on the same sublattice in strained graphene systems. In 
all cases black solid lines and circles (red dashed lines and 
squares) represent the real (imaginary) part of the GF. The 
upper panels represent the GF for a separation of 40 Ia in 
the armchair direction and the lower panels a separation of 
40/z in the zigzag direction. The left (right) panel in each 
case represents the GF for a strain in the armchair (zigzag) 
direction of e = 0.05. An excellent agreement is seen for each. 



the GF, whose corresponding expressions are 
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Q±iE,s) = cos 



\h\it2±E) (4^2 
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2ti 
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(6) 



The strain dependence in these cases enters through ti 
and ^2, given by Eqs. ^ and In Fig [5] we demon- 
strate the remarkable agreement between these expres- 
sions and numerically calculated GFs for a representative 
sample of separation and strain directions. 

RKKY interaction in strained graphene The be- 
haviour of the magnetic coupling can be extracted from 
Eq. dl]) quite easily when the GFs are expressed in the 
form shown in Eq. The integration procedure is 

identical to that for unstrained graphene [39[ and can 
be reduced to a sum over Matsubara frequencies. The 
functions B{E,e) = A^{E,e) and Q{E,e) are expanded 
around Ep and in the low temperature limit T — )■ 0, we 
find 



JsA-lmJ2^^^ cosi2QiEF,e)D) (7) 



where 



{-lYVlB('){E,,s) 
Ji[Ef,s)- — — ZWm — 



is the distance-independent coefficient for the i-th term 
in the power series, ^ is a non-negative integer and 
B^^\EftS) is the £-th order energy derivative of B{E,e) 
evaluated at Ep, resulting from its Taylor expansion. In 
general the leading term in the series should determine 
the decay rate of the coupling. For the undoped case it 
can be shown, for both strained and unstrained graphene, 
using Eqs. © and ^ that the coefficient S(')(0,£) = 0, 
so that the I — 1 term dominates and J{Ep 0) ^ D^^ . 
Thus we should not expect to change the decay rate of 
the interaction by applying uniaxial strain. To study how 
strain does affect the coupling, we examine Eq. ([T]) in un- 
doped graphene as strain is applied and then increased. 
We define the strain-dependent amplification, j3, as the 
ratio between the strained and unstrained couplings, 



Jba{£) 



(9) 



i. Armchair separations The periodicity of the cou- 
pling, determined from Q{E, e) in Eq. ([5]), is clearly inde- 
pendent of ti and t2, and thus strain, for E = 0. Thus the 
only effect that strain can have is a distance-independent 
amplification or suppression arising from the J'l term in 
Eq. ([7]) . In Fig we plot the numerically-calculated 
coupling between moments on the same sublattice as a 
function of armchair separation. This quantity is shown 
for zero strain (black line), and for strains of £ = 0.05 
in the armchair (red, dashed) and zigzag (green, dot- 
dashed) directions. The results, shown in log-log form, 
confirm that the decay rate is unaffected as the strained 
plots lie parallel to the unstrained case. We note that the 
coupling is enhanced by zigzag strain, and suppressed by 
armchair strain. To study the effect of increasing strain 
we calculate the amplification factor /3(e). Using Eq. ([5]) 
we find a simple analytical form for armchair separations 



/3a = 3to 



t2 



^2 



(10) 



i2Q'iEp,e))^+^ 



We note that this expression has the same form in terms 
of ii and t2 for strains in both high symmetry direc- 
tions, but that the ti and t2 values themselves depend 
on the strain direction, as given by Eq. ([2]). In Fig[3jD) 
we plot the analytic expression for /3a as a function of 
strain in both the armchair (red, solid line) and zigzag 
(green, dashed line) strain directions. A monotonic de- 
crease (increase) in the coupling for armchair (zigzag) 
separations consistent with the results in panel a) is ob- 
served. To confirm the analytic prediction, numerical 
calculations of /3^ are performed for a fixed value of sep- 
aration {D = 20 1 a)- Filled and hollow symbols represent 
calculations performed for sites on the same or different 
sublattice(s) respectively. An excellent agreement with 
the analytic predictions is seen in all cases. The iden- 
tical amplification of couplings between same and op- 
posite sublattice sites is explained by the phase factor 



4 




D (y e 

FIG. 3: a) Log-log plot showing the numerical magnetic 
coupling against armchair-direction separation without strain 
(black, solid line) and with e = 0.05 armchair (red, dashed) 
and zigzag (green, dot-dashed) uniaxial strain applied, b) 
Amplification factor Pa as a function of strain in the arm- 
chair (red, solid line) and zigzag (green, dashed) directions. 
Filled and hollow symbols represent numerical calculations 
for same-sublattice and opposite sublattice cases respectively, 
c) Same as a) but for zigzag separation of impurities, d) Am- 
plification factor /3z as a function of armchair strain for same 
(solid line, filled symbols) and opposite (dashed line, hollow 
symbols) sublattice cases. Lines represent the analytic result 
and symbols the numerical calculations, e) Same as d) but 
for zigzag strain. 



between these couplings, which is zero for armchair di- 
rection separations [s^l • For other separation directions a 
more complex behaviour is expected as this phase factor 
is no longer zero, and the phase of the distant-dependent 
oscillations may also be strain dependent. 

a. Zigzag separations In Fig [3]:) we show the mag- 
netic coupling as a function of zigzag separation for the 
unstrained case (black, solid line) and for strains of e = 
0.05 in the armchair (red, dashed) and zigzag (green, dot- 
dashed) directions. A more complicated, non-monotonic 
behaviour than the armchair case is observed. This arises 
due to the Fermi wavevector in the zigzag direction and 
for unstrained graphene the oscillation has a period of 
"ilz- As strain is applied, the wavevector determining the 
oscillation period varies as Q{e) = Q(0) -f 5k{£) where 
5k — cos^^ ("^) ~ amplification factor in the 

zigzag direction is thus 



tg ^Atl - tl cos" {{Q{{)) + 5k{e))D) 
Vsq cos^QiO)D) 



The first part is a distance-independent term similar to 
/3a which gives a monotonic increase (decrease) in the 
coupling for strain applied in the armchair (zigzag) di- 
rection. Viewed with the armchair results, this suggests 
a trend of strain perpendicular (parallel) to the sepa- 
ration direction amplifying (suppressing) the coupling. 
The second part of Eq. ([TT]) accounts for amplification 
due to the change in the Fermi wavevector with strain 
and leads to oscillations in l3z- The analytic expres- 
sion for /3z is plotted in Fig. [3ji) and e) for armchair 
and zigzag strains respectively for both same-sublattice 
(solid lines) and opposite-sublattice (dashed line) cases 
with D = AOlz- The opposite-sublattice results take into 
account the ^ phase shift from the same-sublattice case 
predicted for zigzag separations [s^ . An excellent agree- 
ment is again noted with the numerical calculations rep- 
resented by filled (same-sublattice) and hollow (opposite- 
sublattice) symbols. The oscillations in the coupling, 
which appear as a function of strain, are very interesting 
and may have significant implications for strain-tuning 
of the interaction. Unlike for armchair separations, a 
small difference in the applied strain can tune the cou- 
pling from zero to several multiples of the unstrained 
value. Since the same-sublattice and opposite-sublattice 
couplings are exactly out of phase in this direction one is 
switched off when the other reaches a maximum. For an 
arbitrary non-armchair separation the coupling will have 
characteristic strain values for which one of the couplings 
is zero but the other is not. Thus strain suggests itself 
as a powerful tool, not only to amplify the interaction 
between impurity moments, but also to switch the inter- 
action on and off and to control the interplay between 
impurities on different sublattices. 

Conclusions In this work we have derived analytic ex- 
pressions for the Green function and RKKY interaction 
in graphene for high symmetry directions (armchair and 
zigzag) of both separation and applied uniaxial strain. 
Since GF methods are used to describe a wide range of 
physical properties, our expressions should prove useful 
in the investigation of strained graphene systems. An ex- 
cellent match is found between these analytical expres- 
sions and full numerical calculations. Similarly, the sim- 
ple closed-form expressions describing the amplification 
of the magnetic interaction in a strained graphene sys- 
tem agree with our numerical results. A general trend 
of amplification for strain perpendicular, and supression 
for strain parallel, to moment separation is noted. Also 
noted are oscillations in the amplification as strain is in- 
creased for moments separated in a non-armchair direc- 
tion. This behaviour is again well-captured by our ana- 
lytic approach. Such oscillations suggest the intriguing 
possibility of selectively turning on or off the coupling be- 
tween moments and in particular of controlling the inter- 
and intra-sublattice couplings independently. Since the 
magnetic coupling underpins a wide range of physical 
features, including overall moment formation and mag- 



5 



netotransport response, the ability to fine tune the cou- 
phng with strain may lead to interesting spintronic appli- 
cations. We hope that further investigation of strained 
graphene systems with magnetic impurities will yield a 
diverse range of tuneable properties suitable for device 
application. 
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